Finite  Element  Analysis  and  Experimentation  of  an  Icosahedron  Frame  under 

Compression 

THESIS 

SEPTEMBER  2015 


Mohammed  AlGhofaily,  Captain,  RSAE 

AEIT-ENY-MS-15-S-051 

DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


Wright-Patterson  Air  Force  Base,  Ohio 


DISTRIBUTION  STATEMENT  A. 

APPROVED  EOR  PUBEIC  REEEASE;  DISTRIBUTION  UNLIMITED. 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  refleet  the  official 
policy  or  position  of  the  United  States  Air  Force,  Department  of  Defense,  or  the  United 
States  Government.  This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not 
subject  to  copyright  protection  in  the  United  States. 


AFIT-ENY-MS-15-S-051 


FINITE  EEEMENT  ANALYSIS  AND  EXPERIMENTATION  OF  AN 
ICOSAHEDRON  UNDER  COMPRESSION 

THESIS 

Presented  to  the  Faculty 
Department  of  Aeronautics  and  Astronautics 
Graduate  School  of  Engineering  and  Management 
Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Command 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  in  Aeronautical  Engineering 

Mohammed  AlGhofaily 
Captain,  RSAF 

September  2015 


DISTRIBUTION  STATEMENT  A, 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


AFIT-ENY-MS-15-S-051 


FINITE  EEEMENT  ANALYSIS  AND  EXPERIMENTATION  OF  AN 
ICOSAHEDRON  UNDER  COMPRESSION 


Mohammed  AlGhofaily 
Captain,  RSAF 

Committee  Membership: 

Dr.  Anthony  Palazotto 
Chair 

Dr.  Marina  Ruggles-Wrenn 
Member 


Lt.  Col.  Anthony  M.  Deluca,  PhD 
Member 


AFIT-ENY-MS-15-S-051 


Abstract 


A  nonlinear  analysis  of  a  3D  icosahedron  frame  was  conducted  using  ABAQUS, 
under  a  compressive  load  in  which  collapse  occurs.  The  frame  was  created  in  SolidWorks 
using  the  material  properties  of  the  3D  plastic  building  material  VeroBlue.  Two 
considerations  were  made:  the  computational  features  of  the  frame,  and  the  comparison 
between  the  experimental  solution  to  the  numerical  solution.  Two  studies  were  also 
considered  in  this  research:  The  first  study  was  a  comparison  between  the  seven-inch  and 
the  four-inch  icosahedron  model  with  identical  weight  to  buoyancy  ratios;  and  the  second 
study  was  a  comparison  between  the  seven-inch  icosahedrons  frame,  with  the 
experimental  results  of  the  data  that  was  collected.  The  results  of  both  studies  are 
analyzed  to  compare  against  and  verify  the  results  of  numerical  computational  model. 
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I.  Introduction 


1.1  Chapter  Overview 

The  first  airships,  or  dirigibles,  were  developed  in  the  1800s  [4],  These  simple, 
lighter-than-air  vehicles  (LTAV)  were  filled  with  lifting  gas  to  achieve  positive 
buoyancy  in  air.  Buoyancy  is  the  force  exerted  on  an  object  by  the  fluid  in  which  it  is 
submerged.  The  new  development  of  a  LTAV  which  generates  lift  by  evacuating  the 
air  inside  of  the  structure,  creating  an  inner  vacuum,  had  not  been  researched  at  that 
time.  This  research  introduces  the  development  of  a  numerical  model  for  an 
evacuated  frame  in  the  shape  of  an  icosahedron.  The  frame  supports  the  external 
pressure  acting  on  the  structure,  caused  by  evacuating  the  internal  gas,  which  results 
in  a  vacuum.  The  finite  element  results  are  then  compared  to  the  experiments  of  the 
frame  under  a  compressive  applied  load  (where  the  frame  is  constructed  through  the 
use  of  3D  printing). 

This  chapter  will  describe  the  objective  of  this  research,  and  highlight  the  history 
of  LTAVs,  with  particular  consideration  given  to  why  researchers  decided  upon  the 
icosahedron.  The  assumptions  and  methodology  used  in  this  thesis  will  also  be 
presented. 

1.2  Objective 

The  design  and  manufacture  of  an  icosahedron-shaped  frame  for  a  LTAV,  to  test 
under  compressive  load,  are  of  particular  interest.  The  objectives  of  this  research  are:  to 
gain  a  better  understanding  of  the  compressive  load  resistance  of  an  icosahedron  frame; 
design  and  manufacture  the  frame  using  a  3D  printer;  the  experimental  technique  itself; 
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and  the  evaluation  and  comparison  of  the  finite  elements  analysis  of  the  plastic  material 
used  by  the  3D  model.  An  additional  objective  relates  to  the  weight  to  buoyancy  ratio, 
which  was  found  by  another  researcher  that  it  was  possible  to  establish  a  constant  ratio, 
and  vary  the  size  of  the  structure,  including  the  frame  and  membrane  without  variation  in 
the  internal  stress,  which  was  studied  for  the  frame  itself. 

To  evaluate  the  structure,  a  nonlinear  finite  element  analysis  (FEA)  using  the 
Newton  Raphson  technique  has  been  used.  The  Newton  Raphson  technique  is  a  method 
to  find  successively  better  approximations  of  the  roots  of  a  real  valued  function  [15]. 
Furthermore,  an  eigenvalue-eigenvector  analysis  had  been  carried  out  to  give  a 
comparison  to  the  nonlinear  analysis.  The  Archimedes  principle  states  the  relationship 
which  allows  a  structure  to  become  lighter-than-air  (ETA)  is  when  an  object  is 
submerged  in  any  fluid,  a  buoyant  force  is  expected  equal  to  the  weight  of  the  displaced 
fluid.  The  FEA  will  provide  the  behavior  of  the  response  to  specific  loading  conditions 
which  was  shown  to  be  nonlinear  to  the  ETAV. 


1,3  History  of  LTAV 

An  evacuated  ETAV  is  not  a  new  concept;  it  has  been  used  for  decades.  The 
floating  air  vehicle  dates  back  to  the  1600s  when  an  Italian  monk,  Francesco  Eana  de 
Terzi  proposed  the  "Aerial  Ship  ",  and  stated: 

"The  preceding  inventions  did  not  exhaust  the  ardour  or  the  curiosity  of  the  human 
intellect,  but  have,  rather,  spurred  it  to  seek  how  men  could,  after  the  fashion  of  birds, 
also  fly  in  the  air.  No  one  has;  however,  deemed  it  possible  to  construct  a  vessel  that 
would  travel  in  air  as  if  it  were  supported  on  water.  In  addition,  it  has  not  been  thought 
practicable  to  make  a  machine  lighter  than  the  air  itself,  which  is  necessary  first  to  do  in 
order  to  accomplish  the  desired  end"  [11]. 
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The  main  idea  of  the  statement  is  that  Terzi  wanted  to  invent  objects  lighter  than  air. 

The  hot  air  balloon  was  designed  in  1709,  and  it  was  the  first  LTAV  [8].  The 
concept  involved  heating  the  air  inside  the  balloon  to  reduce  the  density  [3].  This 
principle  led  to  the  creation  of  dirigibles,  in  which  the  gas  used  inside  the  vessel  was 
comprised  of  Helium  or  Hydrogen.  Dirigibles  were  the  first  air  vehicles;  and  there  are 
three  different  types  of  dirigibles:  non-rigid,  semi-rigid,  and  rigid.  The  rigid  and  non- 
rigid  use  a  framework  to  retain  the  vehicle’s  shape  by  use  of  necessary  internal  gas.  The 
semi-rigid  contains  a  partial  framework  to  distribute  the  lifting  loads  [3].  Dirigibles  are 
not  without  their  own  challenges,  however,  including  speed,  control  limitations,  safety, 
and  poor  ground  handling  [10]. 

Currently,  by  using  modem  materials  and  manufacturing  techniques,  a  vacuum 
can  be  created  to  create  positive  buoyancy,  which  results  from  evacuating  the  stmcture. 
The  sphere  is  an  ideal  shape  for  a  vacuum  LTAV  because  it  achieves  the  greatest 
stiffness  with  the  minimum  weight,  but  is  difficult  to  manufacture.  This  research  focuses 
on  an  icosahedron  stmcture,  which  closely  resembles  a  sphere. 


1,4  The  Icosahedron 

The  icosahedron  shown  in  figure  1  is  a  regular  polyhedron  solid  [24].  One  of  the 
advantages  of  the  icosahedron  is  that  it  is  symmetrical,  containing  20  equilateral  triangles 
that  form  the  shape.  The  reason  why  the  bubbles  form  a  shape  of  a  sphere  and  float,  is 
because  it  is  the  ideal  geometric  arrangement  where  the  weight  of  the  bubble  is  less  than 
the  volume  that  surrounds  it;  therefore,  the  external  pressure  causes  lift.  However,  if  we 
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want  to  develop  a  sphere  that  has  an  internal  vaeuum,  we  have  to  inerease  the  thiekness 
of  the  sphere,  thus,  increasing  the  structure’s  weight.  Traditionally,  these  structures  have 
been  designed  with  very  thick  walls  to  resist  buckling.  The  icosahedron  frame  can 
provide  the  symmetry,  spherical  shape,  and  the  necessary  structural  stiffness  and  this 
result  in  uniform  stress  distribution. 


Figure  1 :  Icosahedron  Frame 


1.5  Methodology 

The  FEM  makes  it  easy  to  study  complicated  structures,  which  allow  analysis  of 
structures  like  the  Icosahedron.  The  computer  program  ABAQUS  is  the  finite  element 
tool  used  in  the  analysis  of  the  structure. 
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When  a  compressive  load  is  applied  to  the  structure,  some  design  challenges 
result,  such  as  structural  instability  and  integrity.  Structural  instability  is  the  geometric 
inability  of  the  structure  to  maintain  equilibrium  due  to  the  incapability  to  resist  external 
work  through  internal  energy  [23].  On  the  other  hand,  structural  integrity  is  the  ability  of 
the  structure  to  resist  the  external  load  without  internal  failure.  The  geometric  shape  and 
material  properties  of  the  structure  have  a  relationship  with  both  structural  stability  and 
integrity. 


1.6  Assumptions  and  Limitations 

TEA  includes  assumptions  and  approximations.  Following  assumptions  were  used 
in  this  research; 

1 .  Structural  frame  model: 

a.  The  entire  load  applied  to  the  skin  is  distributed  to  the  frame  with  no  moment. 

b.  Frame  members  act  like  a  beam. 

2.  Manufacturing  of  the  frame 

a.  Minimum  diameter  for  the  beam  is  0. 1 875  inch. 

b.  Maximum  diameter  for  the  frame  is  7  inch. 


1.7  Overview 

•  Chapter  1:  States  the  objectives  of  the  thesis,  the  history  of  FTAVs,  the 
Icosahedron,  and  the  methodology. 
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•  Chapter  2:  Summarizes  the  theory  present  in  the  relevant  literature  related  to  the 
struetural  behavior 

•  Chapter  3:  The  methodology  of  the  development  of  the  different  models  using 
Finite  Elements  analysis  teehniques. 

•  Chapter  4:  Describes  and  provides  the  results  and  the  analysis  of  the  Icosahedron 
models  and  comparison  for  different  material  properties. 

•  Chapter  5;  Provides  the  research’s  conclusion  and  future  work  recommendations. 
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II,  Theory 


2.1  Chapter  Overview 

This  research  will  explore  two  main  areas:  applying  compressive  load  on  the  top 
of  the  icosahedron  structure  considering  the  computational  features  of  the  frame  and 
comparing  a  constant  weight  to  buoyancy  ratio  (W/B)  between  two  frames.  This  will 
evaluate  the  internal  stress  placed  upon  the  two  different  frames,  while  also  comparing 
the  experimental  solution  to  the  numerical  solution.  Previous  research  by  Ruben  Adomo- 
Rodriguez  and  Trent  Metlen  [18]  [22]  established  a  baseline  for  this  research  by 
providing  an  understanding  of  the  static  response  of  the  structure  to  atmosphere  pressure. 
The  purpose  of  this  chapter  is  to  summarize  the  research  which  has  been  carried  out  to 
date  for  the  lighter  than  air  vehicle  (LTAV),  and  the  behavior  of  an  icosahedron  structure. 
The  discussion  starts  with  previous  research  followed  by  nonlinear  analysis,  the  buckling 
and  the  collapse  compressive  failure,  the  manufacture  of  the  frame,  and  concludes  with  a 
summary. 

2.2  Previous  Research 

Ruben  Adomo-Rodriguez  [18]  completed  his  research  concerning  an  icosahedron 
frame  structure,  which  was  originally  developed  by  Trent  T.  Metlen’s  investigation  of 
LTAVs  “to  become  a  viable  method  of  transportation”  [22].  Archimedes’  principle  states 
that  a  body  immersed  in  a  fluid  is  buoyed  up  by  a  force  equal  to  the  weight  of  the 
displaced  fluid  [11].  This  principle  applies  to  both  floating  and  submerged  bodies.  It 
explains  not  only  the  buoyancy  of  ships  and  other  vessels  in  water,  but  also  the  reason 
behind  why  a  balloon  rises  in  the  air,  and  the  apparent  loss  of  weight  of  objects 
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underwater.  In  determining  whether  a  given  body  will  float  in  a  given  fluid,  both  weight 
and  volume  must  be  considered;  that  is,  the  relative  density,  or  weight  per  unit  of  volume, 
of  the  body  compared  to  the  fluid  determines  the  buoyant  force.  If  the  body  is  less  dense 
per  volume  than  the  fluid,  it  will  float,  or  in  the  case  of  a  balloon,  it  will  rise.  If  the  body 
is  more  dense  than  the  fluid,  it  will  sink.  Relative  density  also  determines  the  proportion 
of  a  floating  body  that  will  be  submerged  in  a  fluid.  If  the  body  is  two-thirds  as  dense  as 
the  fluid,  then  two-thirds  of  its  volume  will  be  submerged,  displacing  in  the  process  a 
volume  of  fluid  whose  weight  is  equal  to  the  entire  weight  of  the  body.  In  the  case  of  a 
submerged  body,  the  apparent  weight  of  the  body  is  equal  to  its  weight  in  air,  less  the 
weight  of  an  equal  volume  of  fluid.  Weight  to  Buoyancy  ratio  (W/B)  is  a  concept  that 
establishes  the  buoyancy  in  a  structure  with  respect  to  its  own  weight  [18]. 

In  his  research,  Adorno-Rodriguez  determined  an  equation  for  selecting  a  material 
that  will  satisfy  the  weight-to-buoyancy  ratio  (W/B)  necessary  to  achieve  lift.  His 
calculation  accounted  for  the  atmospheric  effects,  and  is  shown  in  Equation  (1)  [18]. 


W  _  9.5745trVsfem  +  99.098(2c  -  C^yleawJP frame  ^air,i^air,o 


B 


+  ■ 


[2.5362r3  -  ) 

I  nir 


P  ■  T  - 

^  air, air, I 


(1) 


‘  air,o 


Where: 

B  =  buoyancy  of  the  structure 

c  =  beam  thickness-to-radius  ratio  {c=ti,gaml'<^beam) 

Pair.u  Pair.o^  inner  and  outer  air  pressure,  respectively 

R=  air  specific  gas  constant 

r=radius  of  icosahedron  (0.95 1 1  Ibeam) 
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Tair,o>  Tair.i^  inner  and  outer  air  temperature,  respectively 
W=  structure  weight 
Vr=  volume  reduction 

PframeiPskin^  frame  and  skin  densities,  respectively 

He  constructed  three  different  combinations  of  materials  to  plot  the  W/B  for  seven 
different  models.  Figure  2  shows  his  results  for  the  relationship  of  the  applied  pressure  to 
the  max  von  Mises  stresses,  which  is  the  maximum  distortion  energy  criterion  and  is 
often  used  to  estimate  the  yield  of  ductile  materials  [18]  [27].  The  horizontal  lines 
represent  lines  of  positive  buoyancy  indicating  a  threshold,  which  the  applied  stress  on 
the  structure  must  exceed  for  the  structure  to  float  in  air.  The  dashes  lines  represent  the 
yield  strength  of  the  material. 


Figure  2:  Applied  Pressure  versus  Max  Von  Mises  Stress  for  the  Frame  [18] 
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Figure  2  was  conducted  using  a  static  analysis.  The  plot  shows  the  frame  has 
significant  internal  stresses  which  are  above  the  yield  strength  of  the  material,  and  cannot 
achieve  the  required  positive  buoyancy.  However,  in  this  research,  the  material 
considered  is  VeroBlue,  which  is  the  material,  used  in  the  3D  printer  and  the  focus  will 
be  to  manufacture  the  frame  and  compare  the  analysis  with  experimentation.  The  W/B  for 
only  the  frame  is  shown  in  equation  (2). 

fF  ^  ^frameP frame 
^  ^internal  Pair 


Where: 

W/B  =  buoyancy  Ratio 

P Frame’ Pair~  density  of  thc  frame  and  air,  respectively 

^Framer  ^internaf^  Volume  of  the  frame  and  internal  volume,  respectively. 

The  icosahedron  has  equilateral  triangular  faces  and  the  height  of  each  triangle  can  be 
calculated  using  the  Pythagorean  Theorem,  as  follows: 


2 

Figure  3:  Icosahedron  triangle  [28] 
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Where; 


Ibeam^  Length  of  the  beam 


h2  -  i2  _  /beam.2 
'-beam  ^  2 


h2  —  a2 

"•  '■beam) 


V3 

h  —  _  / 

“  2  ^beam 


The  volume  of  the  icosahedron  is  represented  by  20  triangles  of  height  of  r^, 
where  —  OP,  as  shown  in  Figure  4. 


F 

4^ 


A/. 


o  ^ 


•\B 


6o: 


Figure  4;  Equilateral  Triangle  [18] 


Knowing  the  angle  between  sides  is  60°  : 


DP  =  ‘-if^tan30’  =  ^l,„„ 
OC  =  ^cotl8° 

4 


0P2  ^  0^2  _  Qp2  ^  Qp2  _  _  ppy  ^  Qp2  DPy 


(3)  [28] 

(4)  [28] 

(5)  [28] 


(6)  [18] 

(7)  [18] 

(8)  [18] 
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2  _  ileam  Qo  _  ,  V3 


OP^  = 

16 


fy/3  j  V3.  ^  _  /2  .cot^l8°  1 

^beam  g  ^beamj  ~  ^beam\  ^gJ 


ri=OP  =  I 


beam 


jcot^ 

1 

I  16 

48 

—  »  0.7558/ 


beam 


(9)  [18] 

(10)  [18] 


Where  the  area  and  the  internal  volume  is  shown  in  equation  1 1  and  12 


A  = 


2 

beam 


r,  =  20[(i.4): 


/]  “  +  ^^^beam 


(11)  [18] 

(12)  [18] 


The  volume  of  the  frame  with  30  beams  are  shown  in  equation  13 


Vf  =  30nri 


beam^beam 


(13)  [18] 
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Where: 


fbeam^  Radius  of  the  beam 

Then  combining  the  equation  2,  12  and  13  we  get  the  equation  for  the  radius  of  the  beam 


^beam  ^ 


t  g  JPair 


39.07 Alpfy 


(14)  [18] 


Where: 

r=  Radius  of  icosahedron  frame 

2,3  Nonlinear  Analysis 

The  nonlinear  analysis  analyzes  the  behavior  of  the  structure  and  it  is  a  byproduct 
its  stiffness.  The  term  "stiffness"  can  be  used  to  define  the  difference  between  the 
nonlinear  and  linear  structural  behavior.  In  a  linear  structure  analysis,  the  amount  of  force 
per  unit  displacement  is  only  dependent  upon  material  and  geometric  properties  of  the 
structure.  In  a  nonlinear  geometric  analysis,  the  stiffness  is  related  not  only  to  the 
material  and  geometric  properties  of  the  structure  but  also  to  the  structure’s 
displacements  characteristic;  it  can  be  presented  by  equation  (15)  in  which  the  linear 
matrix  is  [/fjyj  and  the  nonlinear  matrix  [/fijJ  which  has  the  displacement  function 
incorporated  [16]. 

[Kij\{Dj]  =  {fj 

[Kij\  =  [Kij\  +  [AKijiu)]  (15) 
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Where: 


K=  stiffness 
F=  foree 

D=  Displaeement 

For  linear  struetures,  if  the  stiffness  is  represented  in  a  simple  form  (k=F/D),  then  for  a 
nonlinearly  responding  strueture,  this  form  must  be  modified.  In  this  research,  we  look  at 
this  type  of  response  [16].  In  geometric  nonlinearly,  the  nonlinearity  is  constantly 
changing.  In  a  linear  scenario  the  stiffness  k  is  equal  to  applied  force  per  unit 
displacement. 

In  a  nonlinear  analysis,  the  stiffness  is  equal  to  the  rate  of  change  of  force  per 
displacement  thus  (k  —  dP/dD).  The  Newton  Raphson  technique  is  the  most  common 
nonlinear  solution  [15],  which  will  be  discussed  in  the  following  section. 


2.2.1  Newton  Raphson  Technique 

The  Newton  Raphson  technique  is  a  method  to  find  successively  better 
approximations  of  the  roots  for  real  value  function  [15].  It  is  used  to  solve  the  nonlinear 
equation.  The  nonlinear  static  equilibrium  equations  is 


KrLq  =  -F{q) 


(16) 


[K  +  iVi(q)  +  N2iq^)^\^q  =  -[K+^N^{q)+  q  +  R 


(17) 


Where: 


14 


K=  linear  stiffness  matrix 


Kj=  tangent  stiffness  matrix 
Ni=  linear  funetion  of  q 
N2=  quadratie  function  of  q 
q=  nodal  displacement  vector 
R=  nodal  loading  vector 

The  Newton  Raphson  algorithm  has  the  iteration  until  the  residual  force  vector  AR 
becomes  sufficiently  small  until  the  equilibrium  equations  are  satisfied.  This  process  is 
repeated  until  it  reaches  convergence.  Figure  5  shows  point  A  and  B,  which  represent  the 
solution  for  load  increment  [16]. 


R 


Figure  5;  Newton  Raphson  Algorithm  [18] 
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2.3  Buckling  and  collapse  failure 


Buckling  and  collapse  are  both  related  to  the  strueture’s  geometrie  design.  The 
difference  between  them  is  collapse  is  charaeterized  by  a  total  nonlinear  analysis.  All 
failure  modes  are  related  to  the  strueture’s  stiffness  [19]. 

2.3.1  Buckling  failure 

Buckling  occurs  when  the  strueture  cannot  recover  from  its  initial  state  of 
equilibrium  after  being  disturbed.  The  state  of  equilibrium  is  when  the  strueture  is  able  to 
move  back  to  equilibrium  [19].  This  is  called  bifureation  loading,  shown  in  Figure  6, 


P 


o 


Figure  6:  Fixed  end  column  [18] 


The  force  applied  to  the  top  end  causes  a  moment  at  point  O,  whieh  will  eause  the 
eolumn  to  bend.  The  stiffness  will  try  to  restore  the  eolumn  to  its  equilibrium  state. 
When  the  load  is  inereased,  it  will  reaeh  a  point  when  the  stiffness  eannot  restore  the 
strueture  to  the  equilibrium  state  whieh  makes  the  strueture  unstable.  At  this  point,  the 
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column  has  buckled,  and  the  load  beeomes  the  eritieal  load.  The  eritieal  bifureation  load 


is  represented  by  equation  (18)  as  shown  in  Figure  6  [19] [2], 

P  =  ^ 

4L2  (18) 


Where: 

I=area  moment  of  inertia 
L=column  Length 
P=  eritieal  load 

A  linear  bifureation  analysis  is  a  proeess  that  estimates  the  eigenvalue,  which 
represent  the  eritieal  (bifureation)  loads.  This  researeh  only  analyze  the  linear  bifureation 
analysis  to  estimate  the  mode  shapes  of  the  strueture,  to  visualize  the  possible  eollapse 
mode’s  shapes.  The  point  when  the  analysis  ereates  the  nonlinear  response,  yielding  a 
maximum  value  without  bifureation,  is  ealled  collapse. 

2.3.2  Collapse  Failure 

Collapse  oeeurs  when  the  strueture  is  not  eapable  of  earrying  its  load.  At  this 
point  the  strueture  is  globally  unstable.  When  the  applied  force  produees  reduetion  in 
stiffness,  and  the  strueture  ean  no  longer  resist  this  load  in  equilibrium,  the  eollapsed 
mode  is  introdueed  [6]. 

The  strueture  eonsidered  beeomes  vulnerable  to  eollapse  when  the  pressure  is 
transferred  from  the  skin  to  the  frame.  Eaeh  member  of  the  frame  is  subjeeted  to  axial 
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load  and  transverse  load.  In  this  researeh,  the  beam  is  a  solid  circular  member  as  shown 


in  Figure  7. 


Figure  7:  Beam  cross-section  for  icosahedron  frame 

2,4  Manufacturing  the  frame 

To  manufacture  the  frame,  the  3D  printer  is  used  as  shown  in  Figure  8.  3D 
printing  or  additive  manufacturing  is  a  process  of  making  three-dimensional  solid  objects 
from  a  digital  file  [25].  The  creation  of  a  3D  printed  object  is  achieved  using  additive 
processes.  In  an  additive  process,  an  object  is  created  by  laying  down  successive  layers  of 
material  until  the  entire  object  is  created.  Each  of  these  layers  can  be  seen  as  a  thinly 
sliced  horizontal  cross-section  of  the  eventual  object  [17].  The  object  is  created  using  a 
virtual  design.  The  virtual  design  is  made  in  a  SolidWorks  Program  in  a  type  of  CAD  file 
(Computer  Aided  Design)  [20]. 

To  prepare  a  digital  file  for  printing,  the  3D  modeling  software  “slices”  the  final 
model  into  hundreds  or  thousands  of  horizontal  layers.  When  the  file  is  uploaded  in  a  3D 
printer,  the  object  can  then  be  created  layer  by  “sliced”  layer.  The  3D  printer  reads  every 
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slice  (or  2D  image)  and  ereates  the  object,  blending  eaeh  layer  with  hardly  any  visible 
sign  of  the  layers,  resulting  in  the  three-dimensional  object  [24]  [21], 


Figure  8;  3D  printer 

There  are  seven  proeesses  to  print  a  3D  model  [26]: 

1.  Vat  Photopolymerisation 

2.  Material  Jetting 

3.  Binder  Jetting 

4.  Material  Extrusion 

5.  Powder  Bed  Fusion 

6.  Sheet  Lamination 

7.  Direeted  Energy  Deposition 
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We  will  only  focus  on  the  material  Jetting. 

Material  Jetting  is  a  process  applied  in  droplets  through  a  small  diameter  nozzle, 
similar  to  the  way  a  common  inkjet  paper  printer  works,  but  it  is  applied  layer-by-layer  to 
a  build  platform  making  a  3D  object  and  then  hardened  by  ultra  violet  (UV)  light  as 
shown  in  Figure  9  [26]. 


Figure  9;  Material- Jetting  schematics  [23] 


With  some  of  the  model  materials,  the  3D  printer  also  jets  a  gel-like  support  material  for 
overhangs  and  complex  geometries.  This  support  material  is  easily  removed  by  hand  and 
with  water  pressure. 
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Figure  10;  Icosahedron  frame  with  support  material 


Figure  1 1 ;  Icosahedron  frame  after  cleaning  the  support  material 
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The  ObjetSOO  ConnexS  3D  printer  is  used  to  print  the  frame.  The  tray  of  this 
printer  is  19.7  x  15.7  x  7.9  in.  It  can  build  multiple  materials  with  multiple  colors  and 
good  surface  finish.  VeroBlue  comes  from  PolyJet’s  family  of  rigid  photopolymers  that 
provide  excellent  detail,  visualization,  durability,  and  strength  [21]. 

In  this  research,  two  studies  of  material  properties  were  considered.  Boston 
University  created  a  dog  bone  and  generated  stress-strain  diagrams  [13].  This  study  is 
very  close  to  the  manufacture  value,  which  is  reported  by  Objet  [17]. The  second  study 
shows  the  values  for  tensile  stiffness  being  lower  than  what  is  reported  by  Objet  [12]. 
Virginia  Tech  investigated  the  Vero White  material  and  the  changes  in  part  tensile 
strength  and  tensile  modulus  due  to  changes  in  orientation  and  part  spacing  [12]. 

Vero  White  has  very  similar  properties  to  VeroBlue. 

Boston  University  did  a  study  of  this  material.  Their  study  is  to  investigate  the 
properties  of  plastics  within  prototyping.  They  created  two-dog  bone  made  of  VeroBlue 
plastic  and  generated  a  Stress-Strain  Diagrams,  representing  the  physical  properties  of  the 
material  [13]. 


Figure  12:  VeroBlue  dog  bone  [13] 
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Figure  13:  Stress-Strain  Diagrams  [13] 


Table  1:  Physieal  properties  of  VeroBlue  [13] 

Elastic  Modulus  Yield  Strength  Ultimate  Tensile  Strength 


2.29  GPa 

50  MPa 

63  MPa 

3.3e5  psi 

7252  psi 

9137  psi 

On  the  other  hand,  Virginia  Teeh  [12]  investigated  the  ehanges  in  a  part’s  tensile 
strength  and  tensile  modulus  due  to  ehanges  in  X-Y  orientation,  Z-orientation,  and  the 
part’s  spaeing  in  the  X-Y  plane,  whieh  is  more  applieable  to  the  experimental  results  of 
this  thesis  due  to  the  spaee  inside  the  ieosahedron.  They  then  analyzed  the  specimens’ 
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tensile  strength  and  tensile  modulus,  using  a  scientific  method  approach,  based  on  three 
parameters:  the  in-build  plane  part  orientation  (X-Y),  the  out-of-build  plane  part 
orientation  (Z),  and  the  distance  between  specimens  [12].  The  results  show  that  the 
spacing  has  the  largest  effects  on  tensile  strength.  This  study  is  used  as  a  guide  to 
compare  to  the  experimental  results  in  this  research.  They  created  an  Ishikawa  Diagram, 
which  is  the  potential  factor  which  could  produce  deviation  in  the  mechanical  properties 
of  the  Vero White  material  as  shown  in  Figure  14  (Cause-and-Effect  /  Fishbone  Diagram) 
[12]. 
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Figure  14:  Ishikawa  Diagram  for  PJD-3DP  Process  [12] 

The  diagram  explains  the  potential  factors  that  can  cause  the  low  and  variable 
tensile  strength  and  tensile  modulus  of  the  Vero  White  material.  Some  causes  include: 
incorrect  geometry,  3D  printing,  methods  used  to  print  the  object,  and  environmental 
differences. 
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Virginia  Tech  analyzed  three  parameters  X-Y  Orientation,  Z  Orientation,  and  Part 


Spacing  as  shown  in  figure  15. 


Parameter:  X-Y  Orientation 


XY  YX 


Description:  The  in-plane  build  orientation  of  the  part  Hypothesis:  Orienting  parts  across  print-head  paths,  may 
such  that  its  length  is  parallel  or  perpendicular  with  lead  to  lower  mechanical  properties  due  to  banding  from 

respect  to  the  front  of  the  build  tray.  discretized  jetting  nozzles. 

Parameter:  Z-Orientation  _ 

YZ  or  XZ  (angled) 

Description:  The  build  direction  orientation  of  the  part  Hypothesis:  An  increase  in  the  number  of  layers  causes 

such  that  the  width,  not  the  length,  lays  flat  or  angled  curing  “print-tlirough”  and  thus  increased  mechanical 

with  respect  to  the  X-Y  plane.  properties  for  parts  with  widths  aligned  in  the  Z  plane 

(as  noted  for  SL  in  [17]). 

Parameter:  Part  Spacmg 


Tight  Spacmg  Far  Spacing 

Description:  The  overall  spacing  between  parts  in  X-Y  Hypothesis:  Smaller  spacing  may  lead  to  increased 
plane  of  the  build-tray.  mechanical  properties  due  to  potential  UV  over-cure 

effect. 


YX  or  XY  (flat) 


Figure  15:  Experimental  Parameters  [12] 

The  researchers  printed  eight  build  trays.  Five  specimens  were  printed  per  build 
tray,  with  each  build  tray  corresponding  to  a  specific  experiment  with  uniform  orientation 
among  all  printed  specimens  [12].  Tensile  tests  were  then  conducted,  with  the  results 
shown  in  Table  2  and  Figure  16  for  mean  tensile  strength  and  modulus  [12]. 
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Table  2:  Mean  Experimental  Data  [12] 


Experiment 

Mass  (g) 

Thickness 

(mm) 

Width 

(mm) 

Tensile 

Strength 

(MPa) 

Tensile 

Modulus 

(MPa) 

YXT 

10.1 

3.23 

13.4 

32.3 

1501 

YXF 

10.1 

3.23 

13.4 

26.0 

1577 

YZT 

10.5 

3.59 

13.0 

31.2 

1696 

YZF 

10.1 

3.43 

13.0 

24.2 

1176 

XYT 

10.0 

3.25 

13.4 

35.3 

1665 

XYF 

10.0 

3.22 

13.3 

29.0 

1719 

XZT 

10.8 

3.66 

13.0 

37.8 

1874 

XZF 

10.4 

3.56 

13.0 

22.9 

1284 

The  first  two  letters  represent  the  orientation,  and  the  third  letter  represent  for  Part 
Spacing  (T  and  F  for  tight  and  far  spacing,  respectively). 
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YXT  YXF  YZT  Y2F  XYT  XYF  XCT  X2F 
Experiments 


Figure  16;  Mean  tensile  Strentgh  (left)  ,  Mean  tensile  modulus  (right)  [12] 
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Considering  the  design  of  the  ieosahedron,  it  has  beams  at  various  angles  and  the 
spaeing  between  these  beams  are  signifieantly  large  in  dimension,  the  most  logieal 
distribution  of  the  modulus  and  the  tensile  strength  will  fall  in  line  with  YZF  or  XZF. 
This  eonelusion  will  be  verified  when  the  eomparison  between  the  FE  model  and  the 
experiment  is  diseuss. 

2.5  Summary 

Nonlinear  analysis  ean  solve  the  eomplieated  mathematieal  problems.  One  of  the 
most  eornmon  methods  to  solve  these  is  the  Newton  Raphson  technique.  There  are  some 
considerations  necessary  when  evaluating  the  structure.  Collapse  is  when  the  structure 
cannot  carry  any  additional  load,  thus,  nonlinear  analysis  will  be  a  main  numerical 
consideration.  In  addition,  the  frame  used  in  this  research  is  manufactured  using  the  3D 
printer.  Also,  two  material  properties  studies  are  being  considered,  with  experimentation 
being  carried  out  to  verily  the  analysis. 
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III.  Methodology 


3.1  Overview 

To  gain  a  better  understanding  of  a  compressive  load  carried  by  an  icosahedron, 
it  is  necessary  to  trace  the  computational  features  of  the  frame,  and  creation  of  the  frame 
model  is  required.  Metlen  and  Adomo-Rodriguez  created  a  frame  model  capable  of 
producing  important  information  about  the  static  response  of  the  icosahedron  [22]  [18]. 
This  chapter  discusses  the  process  used  to  run  the  studies  and  analysis,  model 
development,  and  the  experimental  test  setup  that  is  used  in  this  research. 


3.2  Design  of  the  Frame 

The  reason  for  comparing  the  two  different  size  frames  with  equal  (W/B)  ratio 
was  to  determine  if  the  internal  stress  levels  were  similar  in  each.  Adomo-Rodriguez 
developed  the  icosahedron  model,  and  it  was  the  baseline  for  this  research  [18].  Figure  17 
shows  the  icosahedron  frame  model  used  in  ABAQUS.  Two  icosahedrons  have  been 
designed,  a  seven-inch  diameter,  and  a  four-inch  diameter.  For  an  initial  study  with  the 
goal  of  seeing  if  an  equal  (W/B)  can  produce  a  similar  stress  field  between  frames.  Two 
material  properties  for  VeroBlue  are  used  for  this  research.  One  study  is  to  compare  a 
seven-inch  icosahedron  with  a  four-inch  icosahedron  with  the  same  weight-to-buoyancy 
ratio  (W/B)  [13].  The  second  study  is  to  compare  the  physically  constmcted  seven-inch 
icosahedron  frame  from  a  3D  printed  model  tested,  with  the  FE  model. 

First,  let  us  discuss  the  comparison  of  two  frames  that  have  equal  (W/B)  using  a 
material  utilized  in  the  3D  printer  called  VeroBlue.  The  dimensions,  as  well  as  the  set  of 
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material  properties  for  VeroBlue  studied  in  this  case,  are  listed  in  Table  3  and  4.  The 
finite  element  model  was  created  with  a  weight-to-buoyancy  ratio  of  equal  magnitude, 
utilizing  Equation  (2).  The  radius  of  the  beam  for  the  seven- inch  frame  was  chosen  based 
on  the  minimum  beam  dimension  the  3D  printer  can  process;  less  than  that  may  break  the 
frame  due  to  the  water  jet  used  to  clean  the  support  material  (Figure  18). 
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Water  jet 


Figure  18;  Water  jet  cleaning  off  the  support  material 
Table  3;  Icosahedron  seven- inch  diameter 


Dimension 

Units 

Diameter  (vertex  to  vertex) 

0.1778  (7) 

m  (in) 

Beam  cross  section  Radius 

0.00238125  (0.0937) 

m  (in) 

Length  of  the  beam 

0.09144  (3.6) 

M(in) 

VeroBlue  Density 

1190 

kg/m^ 

VeroBlue  Modulus  of  Elasticy 

2.2945 

GPa 

VeroBlue  Poisson’s  Ratio 

0.35 

Unit  less 
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For  the  four-inch  icosahedron,  the  cross  section  of  the  beam  is  calculated  using 
the  same  W/B  for  the  seven-inch  model.  Using  equation  (2),  (12)  and  (13)  to  calculate  the 
W/B, 


W 

'b 


^frameP frame 
^internal  Pair 


(2) 


Where 

^Framer  ^internaf^  Volume  of  the  frame  and  internal  volume,  respectively. 


Vf  —  307rr^g£j^ 

^beam  (13)  [18] 

Vf  =  307r(0.00238125)(0.09144)  =  4.887xl0-^m^ 

r,  =^(3+V5)(L.„(12) 

r,  =  :2(3  +  V^0.09144=  =  1.668xl0-^m^ 

W  4.887x10-^1190 

—  = - =  28  46 

B  1.668x10-31.225 

Then,  using  this  result  for  the  radius,  the  beam  can  be  formed  using  equation  5 


rbe 


=  r 


t  g  JPair 


39.07  A2pfy. 


^beam 


0.0508  K38-46)(1.225)  ^  Q  QQ;j^3g]^  ^ 
-Vj  39.0742(1190) 


Where  Pair  =  1-255  kg/m^ 


(14)  [18] 
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Table  4:  Icosahedron  four- inch  diameter 


Dimension 

Units 

Diameter  (vertex  to  vertex) 

0.1016(4) 

m  (in) 

Beam  cross  section  Radius 

0.001391  (0.0547) 

m  (in) 

Length  of  the  beam 

0.05232  (2.06) 

m  (in) 

VeroBlue  Density 

1190 

kg/m^ 

VeroBlue  Modulus  of  Elasticy 

2.2945 

GPa 

VeroBlue  Poisson’s  Ratio 

0.35 

Unit  less 

Using  the  Adorno-Rodriguez  finite  element  convergence  study,  the  beams  were 
constructed  with  B32  beam  elements  [18].  B32  beams  in  ABAQUS  are  Timoshenko 
beams  that  allow  for  transverse  shear  deformation,  and  use  a  quadratic  interpolation 
between  nodes  [7].  The  number  of  elements  used  in  this  analysis  was  540  elements. 


3.3  Analysis  Techniques 

The  nonlinear  problem  with  large  displacement  is  to  be  considered  for  the 
analysis.  In  this  research,  the  Newton  Raphson  technique  is  used,  which  was  described  in 
the  previous  chapter.  In  ABAQUS,  the  processor  starts  with  initial  load  increments  and 
the  algorithm  changes  with  the  initial  load  increment  by  25-75%  every  time  for  up  to  five 
iterations  [18]. 
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During  the  process,  three  parameters  are  changed.  These  changes  are  amount  of 
maximum  iterations  per  load  increment  (10000  increments),  the  percentage  of  increment 
change  per  iteration,  and  the  amount  of  equilibrium  calculations  made  before  moving 
onto  another  iteration  (the  initial  increment  size  was  lE-05,  with  a  minimum  of  lE-036 
and  a  maximum  of  0.01).  The  automatic  stabilization,  usually  incorporated  when  a 
membrane  is  part  of  the  structure  is  not  necessary  for  this  analysis.  Eor  further 
comparison,  eigenvalue-eigenvector  analyses  had  been  carried  out. 

3.3.1  Modeling  techniques 

In  this  research,  nonlinear  static  analysis  is  conducted.  The  nonlinear  static 
analysis  is  a  procedure  that  solves  for  the  increment  equilibrium  of  a  structure  given  the 
applied  loads  [18].  The  following  modeling  techniques  and  properties  are  shared  in  all 
icosahedron  models  considered  for  analysis: 

•  Dimensionality:  Two  icosahedrons  models  -  The  first  is  a  seven-inch 
diameter,  with  a  beam  radius  equal  to  0.0937  inches.  The  second  is  a  four- 
inch  diameter  with  a  beam  radius  equal  to  0.001391  inches,  derived  from 
the  same  W/B  of  the  seven-inch  model,  using  equations  (2),  (12),  (13)  and 
(14),  respectively. 

•  Eoad:  The  load,  divided  by  three,  is  applied  at  the  top  of  the  three  vertices, 
as  shown  in  figure  19.  The  load  was  applied  until  it  reached  the  collapse 
mode.  The  seven-inch  model  had  a  load  of  720  N  applied,  and  the  four- 
inch  model  had  a  load  of  255  N  applied  (The  total  load  could  easily  be 
applied  at  a  reference  point  within  ABAQEIS). 
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•  Boundary  Conditions:  Two  boundary  conditions  (BC)  are  considered.  In 
the  first,  at  the  bottom  of  frame,  the  three  vertices  are  fixed  as  shown  in 
Figure  20.  Only  three  of  the  degrees  of  freedom  (DOF)  are  eonstrained: 
U1=U2=U3=0.  The  second  is  at  the  top  where  only  two  DOFs  are 
eonstrained:  U1=U2=0,  allowing  the  vertical  movement  as  shown  in 
Figure  21  (by  the  way,  these  boundary  conditions  were  verified  in  the 
experiment  sinee  the  beams  attaehed  to  the  frames  vertices  rotated). 

•  Mesh:  the  mesh  was  composed  of  B32  beam  elements  for  the  frame. 

•  Analysis:  the  Newton  Raphson  teehnique  was  seleeted  along  with 
asymmetric  matrix  storage  to  improve  the  convergence  history  [7].  A 
linear  buekling  analysis  was  condueted  in  one  of  the  models  in  order  to 
visualize  buckling  mode  shapes. 

Adomo-Rodriguez  generated  the  Python  code  with  a  MATLAB  eode.  The 
MATLAB  eode  was  modified  to  ehange  the  dimension.  All  MATLAB  codes  and 
Python  codes  are  included  in  Appendix  A  and  Appendix  B,  respeetively. 
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Figure  19;  Applied  Load  at  the  top 


35 


Figure  2 1 :  BC  at  the  top 


The  second  study  is  to  compare  an  equivalent  seven-inch  icosahedron  to  the 
experimental  test.  The  3D  printer  model  is  larger  than  the  original  design,  due  to 
imperfection  in  the  printer.  The  new  diameter  is  7.2  inches,  as  shown  in  Figure  22,  and 
the  cross  section  of  the  beam  varied  due  to  the  manufacturing  uncertainty.  In  order  to 
develop  an  average  beam  dimension,  three  diameter  measurements  along  the  length  were 
taken  for  each  beam,  and  the  final  cross  section  diameter  is  the  average  of  the  thirty 
beams  (shown  in  Table  5). 
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Figure  22:  The  new  dimension  for  3D  model 
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Table  5:  Beam  eross  seetion  diameter 


Beam  cross  section  variation  diameter 
(in) 

Average  diameter 
in  (m) 

1 

0.186 

0.188 

0.188 

0.187  (0.004758) 

2 

0.190 

0.188 

0.186 

0.188  (0.004775) 

3 

0.187 

0.187 

0.190 

0.188  (0.004775) 

4 

0.190 

0.189 

0.187 

0.189  (0.004792) 

5 

0.189 

0.189 

0.190 

0.189  (0.004809) 

6 

0.186 

0.187 

0.187 

0.187  (0.004741) 

7 

0.189 

0.186 

0.187 

0.187  (0.004758) 

8 

0.186 

0.186 

0.186 

0.186  (0.004724) 

9 

0.187 

0.186 

0.187 

0.187  (0.004741) 

10 

0.186 

0.187 

0.187 

0.187  (0.004741) 

11 

0.185 

0.185 

0.185 

0.185  (0.004699) 

12 

0.185 

0.186 

0.189 

0.187  (0.004741) 

13 

0.185 

0.185 

0.185 

0.185  (0.004699) 

14 

0.186 

0.187 

0.187 

0.187  (0.004741) 

15 

0.186 

0.186 

0.186 

0.186  (0.004724) 

16 

0.187 

0.187 

0.187 

0.187  (0.004750) 

17 

0.185 

0.185 

0.185 

0.185  (0.004699) 

18 

0.186 

0.186 

0.186 

0.186  (0.004724) 

19 

0.187 

0.187 

0186 

0.187  (0.004741) 

20 

0.185 

0.185 

0.185 

0.185  (0.004699) 

21 

0.187 

0.187 

0.186 

0.187  (0.004741) 

22 

0.189 

0.185 

0.186 

0.187  (0.004741) 

23 

0.187 

0.186 

0.187 

0.187  (0.004741) 

24 

0.189 

0.190 

0.189 

0.189  (0.004809) 

25 

0.186 

0.186 

0.186 

0.186  (0.004724) 

26 

0.187 

0.186 

0.187 

0.187  (0.004741) 

27 

0.190 

0.189 

0.187 

0.189  (0.004792) 

28 

0.187 

0.186 

0.187 

0.187  (0.004741) 

29 

0.185 

0.186 

0.189 

0.187  (0.004741) 

30 

0.187 

0.187 

0.187 

0.187  (0.004750) 

0.187  (0.004742) 
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3,4  Experimental  Test  Setup 

A  compressive  test  was  performed  using  the  MTS  Systems  Corporation  model 
810  Material  Test  System  servo-hydraulic  load  frame,  equipped  with:  15kN  (3.3  kip) 
eapaeity  MTS  model  661.19E-03.  The  MTS  FlexTest  40  digital  controller  with  MTS 
FlexTest  40  Station  Manager  and  Multipurpose  TestWare  (MPT)  application  was  used 
for  the  input  signal  and  data  eolleetion.  Finear  Variable  Differential  Transformer  (FVDT) 
model  FBB-315-PA-100-1,  was  used  to  verify  the  accuracy  of  the  displacement 
measurements  of  the  MTS  maehine.  The  MTS  maehines  have  grips,  not  a  flat  surface  for 
the  placement  of  the  frame.  Thus,  the  manufacturing  of  a  fixture  was  required  in  order  to 
plaee  the  frame  in  fiat  platform.  The  testing  station  is  pictured  in  Figure  23  to  25. 


Figure  23:  MTS  maehine  with  fiat  platform 
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Figure  24:  Icosahedron  frame  with  flat  platform 


Figure  25:  LVDT  attached  to  icosahedron  frame 
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A  high-speed  camera  was  used  to  capture  the  deformation  of  the  beams,  as  shown 
in  Figure  26.  A  BASLER  camera  (model  number  acA2000-165um)  with  a  20-frames  per 
second  (fps)  was  used  to  record  the  beam  deformation. 


Figure  26;  High-speed  camera 


3.4.1  Mechanical  Testing  Procedure 

The  MPT  application  was  opened  within  the  MTS  FlexTest  40  Station  Manager. 
The  procedure  was  run  in  displacement  control,  at  constant  rate  of  0.02  mm/s  (0.047 
inches/min).  Once  the  icosahedron  was  aligned  with  the  fixture,  the  test  was  initiated. 
After  the  test  was  completed,  the  data,  which  was  stored  as  a  .DAT  file,  was  process  in 
Microsoft  Excel. 
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3,5  Summary 

The  purpose  of  this  ehapter  is  to  discuss  the  design  of  the  frame.  The  analysis 
techniques  used  in  this  thesis  are  discussed,  and  afterwards,  the  experimental  test  setup 
with  mechanical  testing  procedure  is  presented. 
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IV.  Analysis  and  Results 


4.1  Overview 

The  modeling  techniques  for  the  analysis  of  an  icosahedron  structure,  and  the 
experimental  test  setup,  were  detailed  in  chapter  three  and  are  repeated  in  this  chapter. 
This  chapter  will  discuss  the  results  of  the  equal  W/B  between  the  seven  and  the  four- 
inch  icosahedrons  with  a  static  applied  load,  and  the  experimental  results  with  a 
comparison  between  the  finite  elements  model  and  an  experiment. 

4.2  Comparison  with  the  same  W/B 

Two  icosahedrons  were  compared  with  the  equal  W/B  between  seven-inch 
diameter  and  four-inch  diameter,  which  made  the  cross  section  of  the  beams  for  the 
seven-inch  greater  than  the  four-inch.  The  load  was  applied  until  the  structure  reached  the 
collapse  mode.  In  this  comparison,  a  linear  bifurcation  analysis  was  conducted  to 
visualize  the  possible  collapse  mode’s  shapes,  and  a  nonlinear  buckling  analysis  was 
conducted  to  visualize  the  buckling.  The  following  sections  discuss  the  linear  and 
nonlinear  analysis. 

4.2.1  Icosahedron  linear  buckling  analysis 

A  linear  bifurcation  analysis  is  a  process  that  estimates  the  eigenvalue,  which 
represent  the  critical  (bifurcation)  loads  and  the  mode  shapes.  This  procedure  is  only  used 
to  estimate  the  mode  shapes  of  the  icosahedral  structure  and  to  visualize  the  possible 
collapse  mode  shapes.  The  FE  model  calculates  many  buckling  modes,  but  only  the  first 
mode  is  relevant,  since  higher  modes  have  no  chance  of  taking  place  because  the 
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structure  will  have  already  collapsed  [27].  The  eigenvalues  were  650  N  for  the  seven-inch 
model  (Figure  27),  and  226  N  for  the  four-inch  model  (Figure  28).  To  carry  out  an 
eigenvalue-eigenvector  analysis  within  ABAQUS,  the  perturbation  method  had  to  be 
used.  In  order  to  do  the  analysis,  the  step  was  changed  to  linear  perturbation,  buckle 
instead  of  (static,  and  general)  to  carry  out  the  evaluation  of  eigenvalue-eigenvector 
analysis. 


U,  Magnitud€ 
r  +1.175«+00 
1  +1.077e+00 
h  4-9.795e-01 
1  +8.815«-01 
h  +7.836e-01 
r  +6.856«-01 
1  +5.877«-01 
1  +4.897C-01 
+3.918«-01 
I  +2.9384-01 
I  +1.9594-01 
+9.7954-02 
>  +0.0004+00 


Step:  Step-2 

Mode  1:  Eigenvalue  =  650.69 
Primary  Var:  U,  Magnitude 

Deformed  Var;  U  Deformation  Scale  Factor;  +1.778e-02 


Figure  27:  Mode  1  in  the  seven- inch  icosahedron 
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Figure  28:  Mode  1  in  the  four-inch  icosahedron 

The  first  eigenvector  for  both  comparisons  between  the  seven-inch  and  the  four- 
inch  eigenvector  were  very  similar.  Even  though  there  was  a  positive  difference  of  250% 
in  the  eigenvalues,  the  sizes  does  not  make  any  difference. 

4.2.2  Icosahedron  nonlinear  buckling  analysis 

A  nonlinear  buckling  analysis  was  conducted,  to  compare  between  equal  W/B  for 
the  seven  and  four- inch  icosahedron,  which  used  the  VeroBlue  material.  The  first  set  of 
material  properties  [13]  of  VeroBlue  were  used  in  comparing  two  icosahedrons  with  the 
same  weight  to  buoyancy  Ratio  (W/B).  The  second  set  of  material  properties  [12]  was 
used  in  comparing  the  experimental  results  with  the  seven-inch  3D  model,  which  was 
larger  than  the  original  design,  as  presented  in  chapter  3. 
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The  load  was  applied  until  the  structure  reached  the  collapse  mode.  For 
the  seven-inch,  the  load  was  720  N  (Figure  29)  and  for  the  four-inch,  the  load  was  255  N 
(Figure  30).  The  load  was  applied  at  the  top  of  the  three  vertices.  It  should  be  noted  that 
the  beam  at  the  loading  and  support  regions  show  bending.  Therefore,  the  boundary 
condition  removing  rotational  resistance  is  correct. 


S,  Mis«s 


Rel.  radius  =  1.0000,  Angle  =  -90.0000 
(Avg:  75%) 


+6.971C+07 

•f6.401e-f07 

•«-5.831e-t-07 

4-5.261e-f07 

+4.691e-t-07 

•f4.122e-i-07 

4>3.552e->-07 

•»^2.982eF07 

•t-2.412e-t-07 

•t-1.842e-f07 

<t>1.272e<t>07 

•(■7.026e+06 

■«-1.328e+06 


Figure  29:  Seven-inch  icosahedron 
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S,  Mises 

Rel.  radius  s  1.0000,  Angle  s  *90.0000 


(Avg:  75<t>b) 

+  6.S01e+07 
-t-6.251e-«-07 
4^5.702e-»-07 
+  5.153e-«-07 
•f4.603e-f07 
■t-4.054e+07 
+  3.504e-f07 
+  2.955e+07 
+  2.406e-4-07 
+  1.856e-t-07 
+  1.307e-t-07 
-l-7.575e-t-06 
+  2.081e-t-06 


Figure  30;  Four-inch  icosahedron 

In  this  comparison,  the  buckle  and  material  failure  were  considered.  The  material 
failed  at  63  MPa  [13].  Figure  31  and  Figure  32  show  the  stress  at  material  failure,  with 
the  reaction  force  at  the  same  failure  stress  for  the  seven-inch  model  and  the  four-inch 
model,  respectively.  The  force  was  the  total  of  the  reaction  forces  acting  at  the  three 
vertices.  Figure  33  and  Figure  34  shows  the  buckling  in  the  icosahedron  for  the  seven- 
inch  and  the  four-inch  frames,  respectively. 
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The  Stress  at  material  failed 


Total  Load  =  885.9  N 


S,  Mises 

Rel.  radius :  1.0000,  Angle :  -90.0000 
(Avg:  75%) 
r  +6.059e+07 
[  +5.554e+07 
[  +5.049e+07 


A. 


Stt(:  )MlBiirSlitic.Ceief  il-ffSUti 
bnainl  BB63i 

Priam  r«RS,<ls*i 

MtrwOfw^U  Mimtai  Suit  (ictiir  - 1 


RF,  Magnitude 
r  +2.953e>02 
[  +2.707e+02 
[  +2.4dle+02 
U2.21SC+02 
[  +1.969e+02 


A 


DrfiraiSlirU  DihnHlMiiScikFiclr-LnDirn 


Figure  3 1 :  Stress  at  material  failed  (left),  Reaetion  foree  at  the  same  stress  (right)  for  the 

seven-inch  icosahedron 


The  Stress  at  material  failed 

S,  Hises 

Rel.  radius  :  1.0000,  Angle  =  -90.0000 
(Avg;  75%) 
r  +5.932e+07 
[  +S.438e+07 
[  ■f4.944e-f07 


Total  Load  =  315  N 


A. 


RF,  Magnitude 
r  +1.0506+02 
\  +9.626e+01 
[  +8.7Sle+01 
[  +7.876e+01 
[  +7.001e+01 


IncriianI  lIZUtpTiBi 
PiinmVirS.4in 

OtfvwiriVprU  O(lmHtiiriSulir«dir*lO0Qr«ni 


Stt*:  \9aiNr'$tadc.Gim<l  nSUti 
iKimal  lIZBiprm:  D.939S 
Priam  Virlir.HiiilMi 
Min>(S*«rU  DtfnialiaiSultFMir'I.n 


Figure  32:  Stress  at  material  failed  (left),  Reaction  force  at  the  same  stress  (right)  for  the 

four-inch  icosahedron 
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The  material  failed  for  the  seven-ineh  model  at  886  N,  and  the  four-ineh  model 
failed  at  315  N,  as  shown  in  Figure  31  and  Figure  32.  The  buckling  for  the  icosahedron 
was  indicated  using  ABAQUS  multiple  plot  states  with  deform  and  undeform  shape,  as 
shown  in  Figure  33  and  Figure  34.  The  beams  started  to  buckle  for  the  seven-inch  and 
four-inch  at  659  N  and  228  N,  respectively.  Figure  35  and  Figure  36  show  the  plots 
between  load  versus  displacement  for  the  seven-inch  and  the  four-inch,  respectively. 


S,  Mises 

Rel.  radius  s  1.0000,  Angle  =  -90.0000 
(Avg:  75%) 


I 


4-1.635e+07 
•fl.503e-f07 
+  1.371e+07 
•t-1.240e-t-07 
•fl.l08e-f07 
+9.764e-f06 
•t-8.447e-f06 
•f7.130e+06 
4^S.814e+06 
•t-4.497e-f06 
+  3.1816+06 
+  1.864C+06 
+  S.476e+0S 


beams  start  to 
buckle  at 
659  N 


Figure  33:  Beams  start  to  buckle  for  the  seven-inch 
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beams  start  to 
buckle  at 
228  N 


S,  Mises 

Rel.  radius  s  I.OOOO,  Angle  s  -90.0000 
(Avg:  75<H>) 


Figure  34;  Beams  start  to  buekle  for  four-inch 


0  OXMCS  OJXn  OjOOIS  0i)02  OtXQS  OiX)3  oms  OtXH  0.004S 

Displacement  (m) 


Figure  35;  Load  vs  Displacement  for  the  seven-inch  model 


50 


0  0.0005  0.001  0.0015  0.002  0.0025 

Displacement  (m) 


Figure  36:  Load  vs  Displacement  for  the  four-inch  model 

In  Figure  35,  the  model  was  linear  until  it  reached  1.2  mm  displacement  then 
started  to  become  nonlinear.  In  Figure  36,  the  model  reached  0.65  mm  until  it  started  to 
become  nonlinear.  The  models  in  ABAQUS  do  not  consider  the  material  failure  criteria, 
which  made  the  structure  stiffer  and  led  to  continuously  rising  load. 

For  the  seven-inch  model,  beams  show  nonlinearity  at  659  N  (148  lb.)  which  was 
close  to  the  eigenvalue  (650  N),  and  the  material  failed  at  881  N  (198  lb.).  The  four-inch 
model  a  beam  start  to  buckle  at  228  (5 1  lb.),  and  the  eigenvalue  was  226  N,  and  the 
material  failed  at  310  N  (69  lb.).  The  study  shows  that  an  equivalent  (W/B)  between 
frame  dimensions  is  not  a  beneficial  way  to  expand  the  design  of  the  icosahedron,  since 
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there  is  virtually  no  similarity  in  stress  between  the  different  frames.  It  was  shown  by 
Cranston  [5]  that  it  is  possible  for  an  entire  icosahedron  with  an  attached  membrane  to 
have  an  equal  stress  distribution,  and  thus  possible  to  design  for  one  (W/B)  ratio  at  any 
elevation. 


4,3  Comparison  FE  model  with  experimental 

Using  Boston  University  material  properties  [13],  the  results  here  were  more  than 
40%  different  compared  to  the  experimental  results  as  shown  in  Figure  37.  The 
difference  in  the  results  led  to  investigate  the  cause.  The  first  step  was  measuring  the 
actual  size  of  3D  model,  which  was  actually  7.2  inch  in  diameter.  The  second  step  was  to 
measure  each  beam  diameter,  which  varied  due  to  the  manufacturing  uncertainty.  The  last 
step  was  to  research  the  material  properties,  which  was  the  second  study  as  detailed  in 
chapter  2. 

In  this  comparison,  the  material  properties  for  the  FE  model  were  varied,  detailed 
in  chapter  2,  Table  2  [12].  The  experimental  test  setup  was  detailed  in  chapter  3,  as 
shown  in  Figure  40  and  Figure  41.  Three  different  elastic  modulus  (E)  were  chosen  in 
this  research:  the  lowest  value  which  was  1.176  GPa,  the  largest  value  which  was  1.874 
GPa  and  the  mid-value  which  was  1.577  GPa.  Eigenvalue  analysis  was  carried  out  to 
estimate  the  buckling  of  the  beams  as  shown  in  Eigures  38  to  40.  Then,  the  displacement 
control  was  applied  to  the  EE  model.  Eor  the  experiment  test,  only  the  load  data  was 
considered.  A  high-speed  camera,  Eigure  42,  was  used  it  to  capture  the  buckling  and 
compare  to  the  EE  model  as  shown  in  Eigures  43  to  45. 
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Figure  37;  Experiment  result  vs  BU  material  properties 


U,  Magnitude 

+  1.175e+00 
+  1.077e+00 
+9.795e-01 

+a.8ise-oi 

•(■7.8366-01 
-(•6.8566-01 
-►5.8776-01 
-(■4.8976-01 
-(-3.9186-01 
■►2.9386-01 
-►1.9596-01 
-►9.7956-02 
-►O.OOOe-^OO 
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Step:  St6p*2 

Mode  1:  Eigenvalue  =  481.98 
Primary  Var:  tl/  Magnitude 

Deformed  Var:  U  Deformation  Scale  Factor:  •►1.8526-02 


Figure  38;  Eigenvalue  for  1.874  GPa 
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Eigenvalue  for  1.577  GPa 
=  405.59 
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Step:  Step-2 

Mode  1:  Eigenvalue  =  405.59 
Primary  Var:  U/  Magnitude 

Deformed  Var:  U  Deformation  Scale  Factor:  +1.852e-02 


Figure  39:  Eigenvalue  for  1.577  GPa 
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Eigenvalue  for  1.176  GPa 
=  302.46 
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Step:  Slep-2 

Mode  1;  Eigenvalue  =  302.46 
Primary  Var:  U,  Magnitude 

Deformed  Var:  U  Deformation  Scale  Factor:  +l.852e-02 


Figure  40:  Eigenvalue  for  1.176  GPa 
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Figure  41;  Icosahedron  Frame  with  flat  platform 


Figure  42:  High-speed  camera 
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Figure  43:  Icosahedron  before  applied  the  load 


Beam  started 
to  buckle 


Figure  44:  Displacement  control  at  2  mm 
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Beams  buckling 


Figure  45:  Displacement  control  at  14  mm 

Figure  46  shows  the  comparison  between  the  FE  model  with  different  material  properties 
and  the  experimental  results.  The  eigenvalues  for  the  three  different  elastic  modulus  (E) 
1.874  GPa,  1.577  GPa  and  1.176  GPa,  were  481.98  N,  405.59  N  and  302.46  N, 
respectively.  These  values  are  the  start  of  nonlinear  behavior  for  each  case  considered. 
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Displacement  (mm) 


Figure  46:  Comparison  between  FE  model  and  experimental 

The  lowest  modulus  is  elose  to  the  experiment  results,  whieh  means  the  material 
properties  of  the  3D  model  is  mueh  lower  than  what  Objet  the  manufaeturer,  reported 
[17].  Comment  related  to  determination  of  this  low  value  of  E  was  made  in  ehapter  2. 


4,4  Summary 

The  eompressive  load  with  the  same  weight  to  buoyancy  ratio  for  the  seven-inch 
diameter  and  the  four-inch  diameter  as  an  equality  function  causes  the  frame  to  get  stiffer 
as  the  dimension  get  larger.  In  addition,  for  an  equal  (W/B)  the  internal  stress  per  frame 
was  variable.  Eor  the  experimental  results,  the  lowest  elastic  modulus  used  in  the  EE 
model  matched  the  experiment  (The  author  carried  out  a  nano  indentation  test  [14]  and 
obtained  an  E  Value  of  1  GPa,  therefore  it  appears  the  lowest  E  value  is  appropriate  in  the 
experiment). 
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V.  Conclusions  and  Recommendations 


5.1  Overview 

The  previous  ehapters  in  this  researeh  have  discussed  the  history  of  LTAVs,  the 
nonlinear  analysis  of  the  icosahedron  frame,  and  the  methodology;  developing  the 
techniques  necessary  for  the  icosahedron  design;  comparison  between  two  equal  weights 
to  buoyancy  ratios,  a  seven-inch  and  a  four-inch  icosahedron;  as  well  as  comparisons 
between  the  seven-inch  icosahedron  with  the  experiment  conducted.  This  chapter  will 
report  the  important  developments  that  transpired  during  the  research,  with 
recommendations  for  future  research. 

5.2  Conclusion  of  Research 

•  The  same  W/B  ratio  made  the  cross  section  of  the  beams  for  the  seven-inch 
greater  than  the  four-inch,  which  led  to  non-similar  distribution  of  stress. 

•  Applying  a  compressive  load  to  the  top  of  the  icosahedron  structure  produced 
local  beam  bending  near  the  loading,  which  removed  certain  rotational  restraints 
for  the  boundary  condition  applied. 

•  Results  showed  that  the  onset  of  nonlinearity  was  very  close  to  the  first 
eigenvalue  of  the  frame. 

•  Experiments  on  the  VeroBlue  material  need  to  be  carried  out  for  all  analyses. 

•  The  dimension  for  the  seven-inch  was  changed  in  the  3D  printer,  and  therefore 
extensive  fine-tuning  on  the  measurements,  not  only  for  the  geometry,  but  for  the 
material  properties,  was  necessary. 


59 


•  The  lowest  modulus  is  elose  to  the  experiment  results,  whieh  means  the  material 
properties  of  the  3D  model  is  much  lower  from  what  the  Objet  manufacturer 
reported. 

•  The  results  from  the  experiments  prove  that  the  method  of  analysis  for  the  frame 
is  correct. 


5,3  Significance  of  Research 

The  nonlinear  static  response  to  an  icosahedron  frame  has  been  evaluated,  and  the 
computer  TEA  model  used  in  researching  the  structure  has  been  verified.  An 
experimental  test  setup  was  developed,  and  a  comparison  to  the  FEA  model  was  made  to 
verify  the  results.  The  nonlinear  behavior  exhibited  by  the  icosahedron  frame  was 
characterized  for  certain  loads  and  boundary  conditions.  This  development  will  help  to 
verify  the  FEA  model  with  the  experimental  test. 


5,4  Recommendations  for  Future  Research 

•  The  comparison  between  two  equally  weight  to  buoyancy  ratio  icosahedrons  were 
used  only  for  the  frame.  Adding  the  membrane  to  the  frame  on  the  icosahedron 
will  result  in  a  more  accurate  response  prediction. 

•  The  experimental  analysis  was  conducted  using  only  the  frame.  Elsing  a  skin 
material  such  as  Kapton  tape  to  replicate  the  membrane,  will  envelope  future 
vacuum  icosahedron  ETAV. 
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•  A  complete  evaluation  of  the  material  properties  produeed  in  any  3D  additive 
manufacturing  process  must  be  earried  out  if  analysis  and  experiments  are  to  be 
properly  compared.  Published  material  properties  values  may  not  apply  to  a 
speeific  situation. 
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Appendix  A:  MATLAB  Codes 


The  MATLAB  codes  are  used  to  generate  the  Python  code  for  icosahedron  model 
in  ABAQUS.  These  codes  were  generated  by  Adorno-Rodriguez  [18].  The  codes  are  to 
setup  the  geometry  of  the  coordinate  icosahedron  frame  and  the  membrane  in  that  frame 
and  the  material  related  to  that  frame.  The  only  change  made  by  the  author  was  in 
(Main.m)  file  to  setup  the  geometry  of  the  vertices  for  the  frame.  This  required  modifying 
the  radius  of  the  icosahedron  (highlight  in  red,  see  page  64)  everything  else  such  as,  the 
geometry  of  the  beams,  the  load  apply  on  the  frame  and  the  material  properties  were 
represented  in  ABAQUS. 

A,1  Seven-inch  icosahedron  model  (Main.m) 


clc;  clear  all ;  close  all 
%  Input 

I.scratch  folder  =  'Temp  Scratch  Files' ; 
I.job.numcores  =  1; 

I.job.memoryusage  =  1024; 
I.job.numGPUs  =  0; 

%  Static  Step  Info 
I. step. buckle  =  0; 

I. step. stabilization  =  1; 

I.step.steptype  =  0; 

I.step.nonlinear  effects  =  'ON' ; 
I.step.increment  method  =  'AUTOMATIC  ; 
I.step.maxnuminc  =  100000000; 

%  Static  Riks 
I.step.initialArclnc  =  0.1; 
I.step.minArclnc  =  le-12; 
I.step.maxArclnc  =  1; 

I.step.maxLPF  =  2; 

%  Static  General 
I.step.initialinc  =  le+00; 

I.step.maxinc  =  1; 

I.step.mininc  =  le-36; 
I.step.stabilizationratio  =  0.05; 
l.step.stabilizationmagn  =  0.0002; 

%  Linear  Buckle 
I.step.bucknumEigen  =  5; 
I.step.buckmaxlter  =  30; 
I.step.bucknumvectors  =  30; 

%  Load 

I.load.dispcontrol  =  0; 
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I. load. d  =  -3e2  ; 

Lload.P=  101325; 

%  Skin  Sections 
1.  section. nostiffnessskin  =  0; 
l.section.membrane  =  1; 

%  Shell  Only 

1. section. skin  section  idealization  =  'NO  IDEALIZATION' ; 
l.section.skin  section  location  =  'MIDDLE  SURFACE' ; 

%  Tie  Constraint 
l.tie.rotations  =  'OFF' ; 

%  Mesh 

l.mesh.skin  element  typel  =  'M3D3' ; 
l.mesh.skin_element_type2  =  'M3D3' ; 
l.mesh.skin  element  shape  =  'TRl' ; 
l.mesh.skin  seed  number  =  18  ; 

2 

l.mesh.frame  element  type  =  'B32' ; 
l.mesh.frame  seed  number  =  18  ; 
l.mesh.rays  element  type  =  'B32' ; 
l.mesh.raysseednumber  =18; 
l.mesh.stiff  element  type  =  'B32' ; 
l.mesh.stiff  seed  number  =  18  ; 

%  Parameters  for  W/B  ratio  calculation 

l.W_B.rho  =  1.2754; 

l.W_B.g  =  9.81; 

l.W_B.skin  =  0.4; 

l.WB. frame  =  0.5; 

l.WB.rays  =  0; 

l.W_B.stiff=0; 

l.W_B.To  =  293.15; 

l.W_B.Ti  =  l.W_B.To; 

l.W_B.Po  =  101325; 

%  Geometry  and  Material  Selection 
%  Material  Selection 
%  rho  nu  E  Sy  ; 

matl  =  [1870  0.3  440e9  3.73e9  ]; 
mat2  =  [1560  0.37  303e9  5.8e9  ]; 
mat3  =  [2700  0.12  757e9  75.7e9  ]; 
mat4  =  [2570  0.33  400e9  3.6e9  ]; 
mat5  =  [1650  0.2  1000e9  10e9  ]; 
mat6  =  [1844  0.18  303e9  0.4e9  ]; 
mat7  =  [2650  0.18  379e9  1.7e9  ]; 
mat8  =  [2800  0.33  738e9  0.14e9  ]; 
mat9  =  [247  0.33  5.76e9  0.024e9]; 
matl0=  [970  0.33  172e9  3.0e9  ]; 

%  Material  Assignment 
matf  =  mat5; 
mats  =  mat5; 
matr  =  mat5; 
matst=  mat5; 

1. materials,  ffamedensity  =  matf(l);  1. materials.  frame_poisson  =  matf(2); 
1. materials.  If  amemodulus  =  matf(3);  1. materials. ffameyield  =  matf(4); 
1. materials. skin  density  =  mats(l);  l.materials.skin_poisson  =  mats(2); 

1. materials. skinmodulus  =  mats(3);  1. materials,  skinyield  =  mats(4); 
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I.materials.rays  density  =  matr(l);  I.materials.rays_poisson  =  matr(2); 

I.materials. rays  modulus  =  matr(3);  I.materials.rays  yield  =  matr(4); 

I.materials.stiff  density  =  matst(l);  I.materials. stiff_poisson  =  matst(2); 

I. materials. stiff  modulus  =  matst(3);  I.materials. stiff  yield  =  matst(4); 

%  Geometry  (icosahedron) 

I.geometry.r  =  0.1778/2;  %  icosahedron  radius,  m;%  This  is  the  only  changes  in  the  code 
I.geometry.rays  =  0; 

I.geometry.  stiff  =  0; 

I.section.hollow_profile_rays  =  1; 

I.section.hollow_profde_stiff=  1 ; 

I.section.hollow_profile  =  1; 

3 

cl  =0.05; 
c2=0.05; 
c3  =0.05; 
for  i  =  1 

I.geometry.  skinthickness  =  I. geometry. r*I.W_B.rho*I.W_B.  skin/... 

(3.77523  *1.  materials,  skindensity) ; 
if  I.section.hollow_profile  ==  1 

I.geometry.  framebeamradius  =  I.geometry.r*sqrt(I.W_B. frame*!. WB.rho/... 
(39.0742*(2*cl-cl''2)*I.materials.frame_density)); 

I.geometry.  framebeamthickness  =  c  1*1. geometry. framebeamradius; 
else 

I. geometry,  framebeamradius  =  sqi't(I.geomefry.r^2*I.W_B.rho*I.W_B. frame/... 
(39.0742*I.materials.frame_density)); 

I.geometry.  framebeamthickness  =  I.geometry.  framebeamradius; 
end 

if  I.section.hollow_profile_rays  ==  1 

I.geometry.raysbeamradius  =  I.geometry.r*sqrt(I.W_B.rays*I.W_B.rho/... 

(1 9. 82*(2*c2-c2^2)*I. materials. rays_density)); 

I.geometry.raysbeamthickness  =  c2*I. geometry. raysbeamradius; 
else 

I.geometry.raysbeamradius  =  sqrt(I.geometry.r''2*I.W_B.rho*I.W_B.rays/... 

(19.82*1.  materials .  ray  sdensity)) ; 

I.geometry.raysbeamthickness  =  I. geometry. raysbeamradius; 
end 

if  I.section.hollow_profile_stiff  ==  1 

I.geometry.  stiff_beam_radius=  I.geometry.r*sqrt(I.W_B.  stiff*!.  WB.rho/... 

(45. 34*(2*c3-c3''2)*I. materials,  stiffdensity)); 

I.geometry.  stiffbeamthickness  =  c3*I.geometry.stiff_beam_radius; 
else 

I.geometry.stiff_beam_radius=  sqrt(I. geometry  .r''2*I.W_B.rho*I.W_B.  stiff/... 

(45.34 *1.  materials,  stiffdensity)) ; 

I.geometry.  stiffbeamthickness  =  I. geometry. stiffbeamradius; 

end 

end 

%  TEA  Analysis 
I. step. buckle  =  0; 

01  =  icosahedron_fea_inner2(I); 
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Appendix  B:  Python  Codes 


Adomo-Rodriguez  [18]  generated  Python  eodes  to  use  in  ABAQUS.  The  only 
modification  was  in  Python  file  (frame.py).  B.l  below  shows  the  deleted  portion  of  the 
code,  which  was  replaced  by  the  Python  file  (Var  icosahedron.py)  as  shown  in  B.2  that 
was  generated  by  the  MATLAB  code  (Appendix  A)  to  generate  the  geometry  of  the 
vertices  for  the  frame.  Similar  modification  was  used  to  analyze  the  other  models. 

B.l  frame.py 


#  Sets  Working  Directory 
os.chdir(path) 

#  Load  variables 
execfile('Var_frame.py') 


B.2  Var  icosahedron.py 


p_array  =  [(0.000000, 0.000000, 0.088900), (0.079515, 0.000000, 0.039757), (0.024571, 0.075623, 0.039757), (- 

0.064329,0.046737,0.039757),(-0.064329,-0.046737,0.039757),(0.024571,- 

0.075623,0.039757),(0.064329,0.046737,-0.039757),(-0.024571,0.075623,-0.039757),(- 

0.079515,0.000000,-0.039757),(-0.024571,-0.075623,-0.039757),(0.064329,-0.046737,- 

0.039757),(0.000000,0.000000,-0.088900),] 


k  = 

[(1,2,3),(1,3,4),(1,4,5),(1,5,6),(1,6,2),(12,7,8),(12,8,9),(12,9,10),(12,10,11),(12,11,7),(7,2,3),(3,7,8),(8,3,4),( 

4,8,9),(9,4,5),(5,9,10),(10,5,6),(6,10,11),(11,2,6),(2,7,11),] 


fc  =  [(0.034695,0.025208,0.056138),(-0.013252,0.040787,0.056138),(-0.042886,0.000000,0.056138),(- 

0.013252,-0.040787,0.056138),(0.034695,-0.025208,0.056138),(0.013252,0.040787,-0.056138),(- 

0.034695,0.025208,-0.056138),(-0.034695,-0.025208,-0.056138),(0.013252,-0.040787,- 

0.056138),(0.042886,-0.000000,-0.056138),(0.056138,0.040787,0.013252),(0.021443,0.065994,- 

0.013252),(-0.021443,0.065994,0.013252),(-0.056138,0.040787,-0.013252),(- 

0.069391,0.000000,0.013252),(-0.056138,-0.040787,-0.013252),(-0.021443,- 

0.065994,0.013252),(0.021443,-0.065994,-0.013252),(0.056138,-0.040787,0.013252),(0.069391,- 

0.000000,-0.013252),] 


np  =  [(0.017348,0.012604,0.072519),(0.057105,0.012604,0.047948),(0.029633,0.050415,0.047948),(- 
0.006626,0.020393,0.072519),(0.005659,0.058205,0.047948),(-0.038791,0.043762,0.047948),(- 
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0.021443,0.000000,0.072519),(-0.053607,0.023369,0.047948),(-0.053607,-0.023369,0.047948),(- 

0.006626,-0.020393,0.072519),(-0.038791,-0.043762,0.047948),(0.005659,- 

0.058205,0.047948),(0.017348,-0.012604,0.072519),(0.029633,-0.050415,0.047948),(0.057105,- 

0.012604,0.047948),(0.006626,0.020393,-0.072519),(0.038791,0.043762,-0.047948),(- 

0.005659,0.058205,-0.047948),(-0.017348,0.012604,-0.072519),(-0.029633,0.050415,-0.047948),(- 

0.057105,0.012604,-0.047948),(-0.017348,-0.012604,-0.072519),(-0.057105,-0.012604,-0.047948),(- 

0.029633,-0.050415,-0.047948),(0.006626,-0.020393,-0.072519),(-0.005659,-0.058205,- 

0.047948),(0.038791,-0.043762,-0.047948),(0.021443,-0.000000,-0.072519),(0.053607,-0.023369,- 

0.047948),(0.053607,0.023369,-0.047948),(0.060233,0.043762,- 

0.013252), (0.067826, 0.020393, 0.026505), (0.040355, 0.058205, 0.026505), (0.023007, 0.070809, 0.013252), (0 

.042886,0.056366,-0.026505),(-0.001564,0.070809,-0.026505),(-0.023007,0.070809,- 

0.013252),(0.001564,0.070809,0.026505),(-0.042886,0.056366,0.026505),(- 

0.060233,0.043762,0.013252),(-0.040355,0.058205,-0.026505),(-0.067826,0.020393,-0.026505),(- 

0.074453,0.000000,-0.013252),(-0.066860,0.023369,0.026505),(-0.066860,-0.023369,0.026505),(- 

0.060233,-0.043762,0.013252),(-0.067826,-0.020393,-0.026505),(-0.040355,-0.058205,-0.026505),(- 

0.023007,-0.070809,-0.013252),(-0.042886,-0.056366,0.026505),(0.001564,- 

0.070809,0.026505),(0.023007,-0.070809,0.013252),(-0.001564,-0.070809,-0.026505),(0.042886,- 

0.056366,-0.026505),(0.060233,-0.043762,-0.013252),(0.067826,-0.020393,0.026505),(0.040355,- 

0.058205,0.026505),(0.074453,-0.000000,0.013252),(0.066860,0.023369,-0.026505),(0.066860,- 

0.023369,-0.026505),] 

pl=(0.000000,0.000000,0.088900) 

p2=(0.079515,0.000000,0.039757) 

p3=(0.024571,0.075623,0.039757) 

p4=(-0.064329,0.046737,0.039757) 

p5=(-0.064329,-0.046737,0.039757) 

p6=(0.024571,-0.075623,0.039757) 

p7=(0.064329,0.046737, -0.039757) 

p8=(-0.024571,0.075623,-0.039757) 

p9=(-0.079515,0.000000,-0.039757) 

pl0=(-0.024571,-0.075623,-0.039757) 

pi  1^(0.064329,-0.046737,-0.039757) 

pl2=(0.000000,0.000000,-0.088900) 

mpl2=(0.039757,0.000000,0.064329) 
mpl3=(0.012286,0.037811,0.064329) 
mp23=(0.052043,0.037811,0.039757) 
mpl3=(0.012286,0.037811,0.064329) 
mpl4=(-0.032164,0.023369,0.064329) 
mp34=(-0.019879,0.061 180,0.039757) 
mpl4=(-0.032164,0.023369,0.064329) 
mpl5=(-0.032164,-0.023369,0.064329) 
mp45=(-0.064329,0.000000,0.039757) 
mpl5=(-0.032164,-0.023369,0.064329) 
mpl6=(0.012286,-0.03781 1,0.064329) 
mp56=(-0.019879,-0.061 180,0.039757) 
mpl6=(0.012286,-0.03781 1,0.064329) 
mpl2=(0.039757,0.000000,0.064329) 
mp62=(0.052043,-0.03781 1,0.039757) 
mpl27=(0.032164,0.023369,-0.064329) 
mpl28=(-0.012286,0.03781 1,-0.064329) 
mp78=(0.019879,0.061 180,-0.039757) 
mpl28=(-0.012286,0.03781 1,-0.064329) 
mpl29=(-0.039757,0.000000,-0.064329) 
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mp89=(-0.052043,0.0378 11,-0.039757) 
mpl29=(-0.039757, 0.000000, -0.064329) 
mpl210=(-0.012286,-0.037811,-0.064329) 
mp910=(-0.052043,-0.03781 1,-0.039757) 
mpl211=(0.032164,-0.023369,-0.064329) 
mpl210=(-0.012286,-0.037811,-0.064329) 
mpl  1 10=(0.019879,-0.061 180,-0.039757) 
mpl27=(0.032164,0.023369,-0.064329) 
mpl211=(0.032164,-0.023369,-0.064329) 
mp711=(0.064329,-0.000000,-0.039757) 
mp211=(0.071922,-0.023369,0.000000) 
mp27=(0.071922,0.023369,0.000000) 
mp711=(0.064329,-0.000000,-0.039757) 
mpll2=(0.071922,-0.023369,0.000000) 
mpl  16=(0.044450,-0.061 180,0.000000) 
mp26=(0.052043,-0.03781 1,0.039757) 
mp610=(-0.000000,-0.075623,0.000000) 
mp61 1=(0.044450,-0.061 180,0.000000) 
mplOl  1=(0.019879,-0.061 180,-0.039757) 
mpl05=(-0.044450,-0.061 180,0.000000) 
mpl06=(-0.000000,-0.075623,0.000000) 
mp56=(-0.019879,-0.061 180,0.039757) 
mp59=(-0.071922,-0.023369,0.000000) 
mp510=(-0.044450,-0.061 180,0.000000) 
mp910=(-0.052043,-0.03781 1,-0.039757) 
mp94=(-0.071922,0.023369,0.000000) 
mp95=(-0.071922,-0.023369,0.000000) 
mp45=(-0.064329,0.000000,0.039757) 
mp48=(-0.044450,0.061 180,0.000000) 
mp49=(-0.071922,0.023369,0.000000) 
mp89=(-0.052043,0.0378 11,-0.039757) 
mp83=(0.000000,0.075623,0.000000) 
mp84=(-0.044450,0.061 180,0.000000) 
mp34=(-0.019879,0.061 180,0.039757) 
mp37=(0.044450,0.061 180,0.000000) 
mp38=(0.000000,0.075623,0.000000) 
mp78=(0.019879,0.061 180,-0.039757) 
mp72=(0.071922,0.023369,0.000000) 
mp73=(0.044450,0.061 180,0.000000) 
mp23=(0.052043,0.037811,0.039757) 

fcl23=(0.034695,0.025208,0.056138) 

fcl34=(-0.013252,0.040787,0.056138) 

fcl45=(-0.042886,0.000000,0.056138) 

fcl56=(-0.013252,-0.040787,0.056138) 

fcl62=(0.034695,-0.025208,0.056138) 

fcl278=(0.013252,0.040787,-0.056138) 

fcl289=(-0.034695,0.025208,-0.056138) 

fcl2910=(-0.034695,-0.025208,-0.056138) 

fcl21011=(0.013252,-0.040787,-0.056138) 

fcl2117=(0.042886,-0.000000,-0.056138) 

fc723=(0.056138,0.040787,0.013252) 

fc378=(0.021443,0.065994,-0.013252) 

fc834=(-0.021443,0.065994,0.013252) 


67 


fc489=(-0.056138, 0.040787, -0.013252) 

fc945=(-0.06939 1 ,0.000000,0.0 13252) 

fc5910=(-0.056138,-0.040787,-0.013252) 

fcl056=(-0.021443,-0.065994,0.013252) 

fc61011=(0.021443,-0.065994,-0.013252) 

fcll26=(0.056138,-0.040787,0.013252) 

fc2711=(0.069391,-0.000000,-0.013252) 

dispcontrol  =  0 
P=101325.000000 
d=-3. 0000006-^02 

hollow_profile  =  1 
framebeamradius  =  8.954073e-04 
framebeamthickness  =  4.477037e-05 
framedensity  =  1650.000000 
frame_poisson  =  0.200000 
framemodulus  =  l.OOOOOOe-i-12 
frameseednumber  =18 
frameelementtype  =  B32 

skinthickness  =  7.280829e-06 
skindensity  =  1650.000000 
skin_poisson  =  0.200000 
skinmodulus  =  l.OOOOOOe-i-12 
skinseednumber  =18 
skinelementtypel  =  M3D3 
skin_element_type2  =  M3D3 
skinelementshape  =  TRI 

skinsectionidealization  =  NOIDEALIZATION 
skinsectionlocation  =  MIDDLESURFACE 
membrane  =  I 
nostiffnessskin  =  0 

raysselect  =  0 
holIow_profile_rays  =  I 
raysbeamradius  =  O.OOOOOOe-l-00 
raysbeamthickness  =  O.OOOOOOe-l-00 
raysdensity  =  1650.000000 
rays_poisson  =  0.200000 
raysmodulus  =  I.000000e-l-I2 
raysseednumber  =18 
rayselementtype  =  B32 

stiffselect  =  0 
hollo  w_profile_stiff=  1 
stiffbeamradius  =  O.OOOOOOe-i-00 
stiffbeamthickness  =  O.OOOOOOe-i-00 
stiffdensity  =  1650.000000 
stiff_poisson  =  0.200000 
stiffmodulus  =  l.OOOOOOe-i-12 
stiffseednumber  =18 
stiffelementtype  =  B32 
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rotations  =  OFF 


modelname  =  'icosahedron-Model' 
jobname  =  'icosahedron-Job' 
jobnameodb  =  'icosahedron-Job. odb' 

#  Step  Information 
buckle  =  0 
stabilization  =  1 
steptype  =  0 
nonlineareffects  =  ON 
incrementmethod  =  AUTOMATIC 
stepname  =  'Nonlinear-Static, General-wStabi' 

#  If  Buckle 
bucknumEigen  =  5 
buckmaxiter  =  30 
bucknumvectors  =  30 

#  If  General 

initialinc  =  l.OOOOOOe+00 
maxinc  =  1  .OOOOOOe+00 
mininc  =  l.OOOOOOe-36 
stabilizationratio  =  0.0500000000 
stabilizationmagn  =  0.0002000000 

#IfRiks 

initialArcInc  =  l.OOOOOOe-01 
minArcInc  =  l.OOOOOOe-12 
maxArcInc  =  1  .OOOOOOe+00 
maxnuminc  =  100000000 
maxLPF  =  2.000000e+00 

numcores  =  1 
memoryusage  =  1024 
num  GPUs  =  0 
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